Effect of Composition Changes on the Structural Relaxation of a Binary Mixture 
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Within the mode-coupling theory for idealized glass transitions, we study the evolution of struc- 
tural relaxation in binary mixtures of hard spheres with size ratios 8 of the two components varying 
between 0.5 and 1.0. We find two scenarios for the glassy dynamics. For small size disparity, the 
mixing yields a slight extension of the glass regime. For larger size disparity, a plasticization effect 
is obtained, leading to a stabilization of the liquid due to mixing. For all 5, a decrease of the elastic 
moduli at the transition due to mixing is predicted. A stiffening of the glass structure is found as 
is reflected by the increase of the Debye- Waller factors at the transition points. The critical am- 
plitudes for density fluctuations at small and intermediate wave vectors decrease upon mixing, and 
thus the universal formulas for the relaxation near the plateau values describe a slowing down of 
the dynamics upon mixing for the first step of the two-step relaxation scenario. The results explain 
the qualitative features of mixing effects reported by Williams and van Megen [Phys. Rev. E 64, 
041502 (2001)] for dynamical light-scattering measurements on binary mixtures of hard-sphere-like 
colloids with size ratio 5 = 0.6. 

PACS numbers: 64.70.Pf, 82.70.Dd 



I. INTRODUCTION 



The study of glass-transition phenomena in dense col- 
loidal suspensions has received much attention during the 
past years. Such systems are well suited for a test of theo- 
ries since the particles' properties can be tuned within a 
broad range. In particular, one can produce mixtures 
with particles of different sizes and observe effects of 
changing composition or size ratio. Recent experimen- 
tal work of Williams and van Megen jj| has shown that 
even in the simplest such systems, namely binary hard- 
sphere mixtures (HSM), interesting mixing phenomena 
appear for the dynamics close to the glass transition. 
Three effects have been reported going over from a one- 
component to a binary system containing up to 20% by 
volume of smaller spheres: (i) a shift of the glass tran- 
sition to higher packing fractions, (ii) an increase in the 
plateau values of the correlation functions at interme- 
diate times, connected to an increase in the glass form 
factors, and (iii) a slowing down of the initial part of the 
relaxation towards this plateau. 

In this paper, mixing effects in binary HSM are inves- 
tigated in the framework of the mode-coupling theory of 
the idealized glass transition (MCT). The study of glass- 
transition phenomena in colloidal suspensions that are 
good realizations of one-component hard-sphere systems 
has revealed that MCT describes much of the experimen- 
tal facts in these cases || . MCT makes general predic- 
tions for all glass-forming systems, independent of their 
underlying microscopic properties, be they one- or multi- 
component systems. Thus, a universal glass-transition 
scenario has been established, involving scaling laws and 
power-law variations of time scales. These properties 
have been found in many, not only colloidal, systems, as 
reviewed in Refs. j|, But MCT is also able to derive 
detailed results depending on the specific interactions of 
a system. The aforementioned hard-sphere colloids are 



a paradigmatic example for which, among other things, 
the wave-vector dependence of the Debye- Waller factors 
in the glass state has been evaluated and compared with 
results from scattering experiments |B| . The quantitative 
study of model systems allows one to predict general, 
while nonuniversal, trends that arise in certain classes 
of glass formers. Such project has been carried out for 
molecular liquids, where the known differences of reori- 
cntational relaxation for angular momentum I = 1 and 
i = 2 could be explained || The work presented 
here in a similar way aims to explain the general trends 
occurring in a mixture when changing its composition 
or the size disparity of its constituents. Our discussion, 
motivated by the cited light-scattering experiments E), 
focuses on binary HSM with not too large size disparity 
in the species, close to the glass-transition density. 

For a derivation of the MCT for mixtures, the reader 
is referred to Ref. ||. The theory has already been ap- 
plied to analyze conrputer-simulation data for a binary 
soft-sphere mixture [|9|, a binary Lennard- Jones mixture 
Jio| , [lH , a molecular-dynamics model of a silica melt |L2| , 
and of a two-component metallic melt [fl3|| . Also, prop- 
erties of binary HSM in the limit of large size disparity 
|jLH p"5[ |l6| , [l7| , |l8f and of charged hard spheres, par- 
ticularly in their low-density region JO], [2(], have 
been studied in the framework of MCT. Mixing effects 
in a binary HSM have been addressed recently using a 
standard liquid-state mode-coupling approximation, al- 
beit for states of such low density that glassy dynamics 
does not occur [^2[ . MCT equations for mixtures have 
been derived recently within a nonlinear-hydrodynamics 
theory p4| . The found equations are very different from 
the ones analyzed here, and a connection of their implica- 
tions with the light-scattering data B was not discussed. 

The paper is organized as follows. In Sec. |l|, we sum- 
marize the basic formulas specifying the model under 
study. Sections III and |y| discuss our results for the 
fluid-glass transition diagram and the glass-form factors, 
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respectively. We demonstrate in Sec. [v] that these lead 
to two qualitatively different scenarios for the dynamics 
close to the glass transition. Section |v| summarizes the 
results. 



II. DEFINITION OF THE MODEL 

A. General Equations of Motion 

A classical S'-component fluid of N spherical par- 
ticles shall be considered. The fluctuations of the 
partial number-densities shall be denoted as g a (q) = 
J2k exp[igfk ]/\/iV, a — 1, 2, . . . , S, where the sum runs 
over all N a particle positions belonging to species 
a. From this, the partial density correlators are con- 



structed, written <fr a p(q,t) — (g a (q)\gp(q,t)). Here, 
(A\B) = {5 A* SB) with 5 A = A - (A) denotes a scalar 
product in the space of dynamical variables. Angle 
brackets indicate canonical averaging for temperature 
T. The time evolution is generated by a Liouvillian C: 
g a (q,t) = exjp[i£t]g a (q). Since the g a are spatial Fourier 
transforms of a real density variable in an isotropic, trans- 
lational invariant system, $> a p(q,t) is real, even in t, and 
it depends on the wave vector only through q — \q\. An 
evaluation of the density correlators & a p (q, t) is the ma- 
jor aim of this paper. 

The starting point of the theory is the exact reformu- 
lation of the equations of motion using the Zwanzig-Mori 
technique. Considering the limit of a system of colloidal 
particles undergoing Brownian dynamics, this equation 
reads 



It is to be understood as a matrix equation as is 
indicated by the bold symbols. S(q) is the ma- 
trix of partial structure factors defined by S a p{q) — 
(g a (q)\gp(<t)) ■ The short-time behavior of the corre- 
lators is given by = S(q) - r(q)- 1 \t\ + 0(t 2 ), 
where T(q) is a symmetric positive definite matrix 
of relaxation times. It shall be specified in terms 
of short-time diffusion coefficients D° as r a p(q) = 
l/(q 2 D^)5 a p- The memory kernel M(q,t) is given 
through the so-called fluctuating forces, M a p(q,t) = 
{k B T) 2 (x a /m a )(xp/mp)( Q£j Q (q)\K'{t) Q£jp (q)). Here, 
Q is the projector perpendicular to the number densities, 
Q a (q), and the longitudinal parts of the number current 
densities, qj a (q) = Cg a (q). 72.'(i) = exp[iQCQt] is the 
reduced evolution operator. The m a are the masses of 
the species labeled by a, and x a = N a /N are the number 
concentrations. 



Equation (la) is complemented by an approximate ex- 
pression for the memory kernel. The MCT approxima- 
tion for this quantity follows from a straightforward gen- 
eralization of the one-component case ]8(] and gives the 
polar form 



M(q,t)=F[&(t),$(t)}(q) 



(lb) 



r( g )*( 9 ,t) + S( (Z )- 1 *((7,i) 



M(q,t-t')4>(q,t')dt' = 0. (la) 



of a symmetric bilinear form of the density correlators 



^[*«,$( 2 )](9) 



2q 2 x a xg ^— ' ^ 



E E v °«°" (£ w^V' (p)Vw0" % Z p) 



(lc) 



r 



Here, p is the total number density, p = q — k, and 
V a /3f{q,k,p) are vertices quantifying the coupling of a 
force fluctuation of wave vector q to density-fluctuation 
pairs with wave vectors fc and p, respectively. The ver- 
tices are given by the equilibrium structure of the system 
in terms of the Ornstein-Zernike direct correlation func- 
tion c a p(q) and static three-particle correlations. The 
latter shall be expressed in terms of the (g) using the 
convolution approximation. One thus arrives at 



>(q,k,f) = (qk/q)c aa '(k)S aa » + (qp/q)c aa »{p)6 



(Id) 



A few remarks to these equations might be in order. 
The solution to Eqs. ([[]) exists for all t > and it is 
uniquely determined by the initial conditions. For a sys- 
tem with colloidal short-time dynamics, the correlation 
functions <&(q, t) are completely monotone functions p5[ . 
This property is preserved by the specified MCT approx- 
imation. In detail, it implies the following. The matrices 
&(q,t) are positive definite, written <&(q,i) >z 0, for all 
times t and at all q; and for the time derivatives, there 

holds (— l) l d\&(q, t) >z for all I = 1, 2, Furthermore, 

the solution depends smoothly on S a p(q), T a p(q) and 
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F 

-l 



= S(q) 
= 1,2,..., with 



Vap-y^q, k,p) for any fixed finite time interval. At t — > oo, 
bifurcations in the long-time limit F(q) = lim^,^ <f>(q, t) 
may occur at critical points called glass-transition singu- 
larities. The F a p(q) are called the glass form factors and 
they are solutions of the equation 

F( q ) = S(q)-[S( q )- x +F[F,F]( q )]- 1 . (2a) 

In particular, the correct F(q) can be determined 
through an iteration scheme. 

Siq)- 1 +T[F {n - 1 \F^- 1} ](q) 

starting point F (0) (q) = S(q). The sequence F {n) (q) 
converges monotonically towards F(g). For sufficiently 
small vertices, one has the liquid solution F(q) = 0, 
while the glass is characterized by F(q) ^ 0. Solutions 
corresponding to critical points shall be denoted by a su- 
perscript c. To understand the bifurcation scenario, one 
needs to discuss the critical eigenvector of the lineariza- 
tion of Eq. (pa]), H(q), given through 

H(q)-2(S c (q)-F c (q))T c [F c ,H}(q)(S c (q)-F c (q)) = 0. 

This eigenvector is nondegenerate, which implies that all 
MCT bifurcations belong to the type Ag, i — 2, 3, . . ., 
introduced by Arnol'd (g6|. Also, one can choose H(q) >~ 
0. All of the preceding remarks apply to the general case 
of matrix-valued MCT equations; they are not affected by 
the precise form of the vertices entering Eq. (|k). Proofs 
of the cited mathematical properties of Eqs. ([[) and (||) 
can be found in Ref. . In this paper only the simplest 
bifurcations, i.e. those of type A^ are discussed, where a 
jump occurs in F(q) from to the critical value F c (q) y 
0. 

There is an important implication of the cited general 
results that can be substantiated in complete ana logy 
to the one obtained for one-component systems |H, ^9[ . 
There exists a time scale to such that one can express the 
solution of Eq. (|l|) for t ^> to in the form 

*(q,t) = S 1 '*(q)<Kq,i)S 1 ' 2 {q), (3a) 

where t = t/to- The time scale to depends smoothly on 
S a p{q), T a p(q), and V a 0~ f (q,k,p). The sensitive depen- 
dence of <fr on the control parameters is described by the 
completely monotone function q>(q,t). The latter is de- 
termined, up to some arbitrary time scale £», through the 
equation 



(j)(q,t) = m{q,t) 



d_ 

dt jo 

obeying the initial condition 

lim(i/**) 1/3 <K<?>*) = 1. 



m(q,t- t')(f>(q,t')dt' , (3b) 



(3c) 



Here the kernel m(q,t) is given by the mode-coupling 
functional, Eq. (fTc]), as 

m{q,t) = S 1 ' 2 {q)x 
x T [S^VWS 1 / 2 , S 1 ' 2 ct>{t)S 1 ' 2 \ (q)S 1/2 (q) . (3d) 



This means that <j>(q, t) is determined by the quantities 
entering T only, i.e. by the parameters specifying the 
equilibrium structure. In this manner, MCT justifies the 
concept of structural relaxation as opposed to, e.g., tran- 
sient relaxation. The latter, be it Brownian or Newto- 
nian, merely enters the scale to- In particular this implies 
that details of T a p (q) do not affect the structural relax- 
ation apart from influencing to- 

Let us also note the formulas for the longitudinal and 
transversal elastic moduli of the mixture. They are given 
through Green-Kubo relations involving the total mass 
currents ||(], One defines the projector Qhd as 

projecting out g a (q) and the longitudinal mass current, 
= J2 a m aja(q), together with the corresponding re- 
duced resolvent TZ^^z). The latter is the Laplace trans- 
form for frequency z of the corresponding evolution op- 
erator 7^y D (t) = expfiQHD^QHD^l- For the longitudinal 
viscosity, this yields |30j 



V 



L = lim lim 



1 



^09^0 fc B T q 2 



(QHDCJ(q)\TZ' HD (z)Q HD CJ(q)) 



(4) 

At the bifurcation singularity, a nontrivial long-time 
limit F c (q) implies a pole at zero frequency, & c (q, z) ~ 
—F c (q)/z. According to this, the longitudinal modu- 
lus shows a discontinuity 5M£ at the glass transition. 
Rewriting Eq. (||) in terms of the MCT projector Q, one 
gets 

SMI = (gk B T) lim V x a ^[F c , F c ]{q)x p . (5) 

a/3 

The glass, different from the liquid, is characterized by a 
finite shear modulus. A formula similar to Eq. (|^) is ob- 
tained for the shear modulus at the glass transition, 



M% = (gk B T) lim V x a &[F°, F c ](q)x . (6) 

a— »0 — 4 y 

a/3 

Here, the MCT expression for the transverse fluctuating- 
force kernel ^" T (g) is obtained from ^(q) by replacing in 
Eq. (|lc]) the vertices by 

HL'q"(9> k,p) = (q T k/q)c aa/ (k)S aa » 

+ {q T P/ q)Caa" (P) 5 aa' ■ (7) 

In this formula, q 1 - is a vector of length q perpendicular 
to q. 



B. The Binary Hard-Sphere Mixture 

The general theory shall be applied to binary hard- 
sphere mixtures (HSM), consisting of large (A) and 
small (B) particles. If d ai a = A,B, denote the parti- 
cle diameters, the packing fractions of the species read 
ip a = (ir/6)(x a g)da, and the total packing fraction is 
given by <p = <pa + <Pb- The thermodynamic state is 
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characterized by three control parameters. Let us choose 
them to be the total packing fraction ip, the size ra- 
tio S — ds/dA < 1, and the packing contribution of 
the smaller species xb = V'b/V- Whenever composition 
changes are considered in the following, a variation of xb 
for fixed if and d is to be understood. This in turn im- 
plies the number concentration of the small particles to 
vary as 



xb 



x B /S 3 



l + x B (l/S 3 -l) 



(8) 



The procedure is somewhat in between a true addition, 
which would increase both total density and total packing 
fraction, and a replacement of a certain amount of large 
spheres by the same amount of smaller ones, which would 
reduce the total packing fraction. For sufficiently small 
(5, there appears a percolation threshold for the motion of 
the small particles in the glass formed by the large ones. 
This transition and its precursor phenomena shall not be 
considered in this paper. 

Static structure input for our model is taken from the 
Percus-Yevick (PY) approximation ||^, [53). More accu- 
rate solutions of the Ornstein-Zernikc integral equations 
for hard-sphere mixtures are available. Yet, one knows 
from the one-component MCT that improvements aim- 
ing at, for example, thermodynamic consistency have lit- 
tle influence on the glassy dynamics. Unfortunately, the 
quality of the PY approximation at the desired high pack- 
ing fractions is unknown. It is known that large errors 
of the structure factor can occur if one goes over to large 
values of 1 /S [M , but this case is excluded from our dis- 
cussion. 

With the structure factor and the direct correlation 
functions given, the vertices in Eqs. (Id) and (Q) are well 
defined functions of the wave vectors and matrix indices. 
Hence, also the mode-coupling functional T in Eq. ( |Tc|) is 
defined as a triple integral over the components of k; and 
the same holds for the functional T T . After introduction 
of bipolar coordinates and using rotational symmetry, the 
k integrals are transformed to double integrals over k = 
\k\ and p = \q — k\. After performing the q — > limit 
in the functionals, the zero-wave-vector limits entering 
Eqs. (§) and (§) are reduced to one-dimensional integrals 
over k. As a next step, the wave vectors are reduced 
to points on a grid of M values. The grid is chosen as 
qdA = qo + qAq, with qo = 0.2, Aq = 0.4, and q — 
0, 1, . . . , M — 1, unless otherwise stated. The integrals 
are replaced by Riemann sums. The resulting formulas 
are the same as explained explicitly before for the one- 
component systems ]3^] , but additional sums over matrix 
indices occur. As a result, the cited equations refer to 
ones for sets of M matrix correlators, where q serves as 
a label for the correlators. To complete the specification 
of the equations, the short-time diffusion constants are 
taken according to Stokes' law, -D° = C/d a . The unit of 
time is chosen so that C = 0.01. 

For the simple hard-sphere system, it was found that 
a choice of M — 100 is sufficient to avoid cut-off and 



discretization effects for the results [ p5| . Our work was 
done mostly with M — 200, implying a cut-off wave vec- 
tor q*dA = 79.8. This enables us to handle size ratios 
5 > 0.5 with the accuracy used earlier for the simple 
system. Thus, Eqs. (|h^-(|Tcj) are 600 coupled integro- 
differential equations for 600 correlators, and Eq. ( J2a| ) 
formulate 600 implicit equations for the 600 glass form 
factors. The latter equations are solved by the iteration 
mentioned above. Thereby, one gets the form factors, 
discussed in Sec. |y|. Shifting the value for the packing 
fraction tp, one identifies the glass transition points. The 
linear equation for the nondegenerate eigenvector H(q), 
Eq. (^), is solved by a standard routine. This eigenvec- 
tor is us ed to calculate the asymptotic solutions discussed 
in Sec. VB . The critical glass form factors are substi- 
tuted into Eqs. (^) and (|^) so that the contributions to 
the moduli discussed in Se c. |IIl| can be calculated. 

The closed set of Eqs. (|la])-(px[) as well as Eqs.(|3a|) 
([k]) for the time dependence of the correlators are solved 
by a method adapted to this special kind of Volterra 
problem. The solutions of Eqs. ([[]) yield the correlators 
to be compared with the cited light scat terin g data. The 
solutions of Eqs. ( pb| ) are used in Sec. VB to separate 
structural relaxation from the transient dynamics. To 
proceed, one introduces a grid on the time axis of equal 
step size h consisting of N points. The time derivatives 
in Eqs. ( |la| ) and (3b) are replaced by difference relations 
and the convolution integrals by Riemann sums. One 
gets a recursion relation determining the solutions for 
time t n+ i = h(n + 1) from the values ti, I < n. The 
initial values are taken from the short-time asymptote, 
given by <&(q,t) — S(q) ~ — r(q)~ l t or Eq. (|3c[), respec- 
tively. Because of the scale invariance of the equations 
of structural relaxation, the value of does not mat- 
ter. One has to make sure that the results obtained for 
< t < ijv remain stable within the desired accuracy 
upon doubling h or halving N . Then one carries out a 
decimation by setting h i— ► 2h. Thereby, one solves the 
equations up to It n . This procedure is repeated until the 
correlators reach their long-time asymptote. Details are 
explained, e.g., in Ref. ]36| . Typically, the figures in this 
paper have been calculated with h = 10 -6 and N = 256. 
Let us point out that Eqs. (|l|)~(0) are completely analo- 
gous to the ones discussed repeatedly in previous studies 
of one-component systems. Also the numerical meth- 
ods applied for the solution are the ones used earlier for 
the simple case. The additional complication here is the 
handling of matrices; but this is straightforward though 
numerically more demanding. 



III. TRANSITION DIAGRAM 

Cuts through the liquid-glass transition surface in the 
three-dimensional control-parameter space for the binary 
HSM are depicted in Fig. I]. To assure that the results do 
not seriously depend on the discretization used, we show 
as well the glass transition points calculated for S = 0.6 
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FIG. 1: Liquid-glass transition diagram of a binary hard 
sphere mixture (HSM) for size ratios 8 — 0.5 (triangles), 
8 = 0.6 (diamonds), 8 = 0.7 (squares), and 8 = 0.8 (circles), 
plotted as critical total packing fraction tp c versus packing 
contribution of the smaller species, xb ~ (Pb/<P- Full lines are 
guides to the eyes. The dashed lines indicate results obtained 
from tripling the number M of grid points from M — 200 
to M = 600, and the dotted line for 8 = 0.6 shows results 
obtained using M = 100; see text for details. 



and 0.8 from the model with M = 600, Aq = 0.4/3, 
go = 0.2/3. In addition, for 6 = 0.6 the dotted line 
exhibits the result calculated with cutoff q*d,A = 39.8 
and M = 100 which are the discretization parameters 
used in Ref. jj5|. One infers that for xb < 0.3, this 
discretization would be sufficient to produce reasonable 
results. 

For fixed size ratio S < 0.65, the critical packing frac- 
tion first increases upon increasing xb- Since xb = 
and Xb = 1 both represent monodispcrse hard-sphere 
systems, one gets ip c (xB = 0) = <p c {£b = !)• Thus, 
the liquid-glass transition lines for S < 0.65 exhibit a 
maximum at some intermediate values of xb- It is well 
understood that big particles moving in a liquid of much 
smaller ones experience an effective attraction |37j that 
is of purely entropic origin. Such a short-ranged attrac- 
tion leads to a stabilization of the liquid phase, as was 
explained earlier [^8[ |39| |40) . Our result is a direct anal- 
ogon of this depletion- attraction effect. Similarly, from 
the discussion of polymer melts it is known that the intro- 
duction of smaller components into the system typically 
decreases the viscosity, i.e., drives the system further into 
the liquid phase; an effect sometimes called "plasticiz- 
ing". Therefore, the effect found here is an entropically 
induced plasticization effect. 

For less-disparate-sized mixtures, the theory predicts 
an inversion of the effect described above. An exam- 
ple is shown in Fig. [I] for 6 = 0.8, where a decrease of 
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FIG. 2: Maximum relative increase and decrease of the criti- 
cal packing fraction, Aip c (S), according to Eq. (^|), as a func- 
tion of the size ratio S (crosses), together with experimental 
data for random loose packing (triangles, reproduced from 
Ref. cf. text). For the MCT critical packing fraction val- 
ues, two symbols are noted for those S, where a maximum and 
a minimum different from the 8 — 1 value could be identified. 



ip c with increasing xb up to some minimum point is ob- 
served. This is in accordance with similar MCT results 
for a binary soft-sphere mixture ||. It means that the 
introduction of disorder due to a small polydispersity of 
the particles stabilizes the glass state. The transition dia- 
gram is not symmetric with respect to xb — * (1 — xb); our 
theory predicts for 0.65 < S < 0.8 "S"-shaped transition 
lines. 

To get another view on the transition diagram, let us 
define the relative change of tp c with respect to the one- 
component case through 



a^(5) = (^(5,4)-^o)M. 



(9) 



Here, x B are the points at which a maximum or a min- 
imum occurs in ^ c (ib) for fixed 5; <p§ w 0.5159 is the 
critical packing fraction of the one-component system. 
The resulting values are plotted in Fig. g together with 
data taken from Ref. There, results for Aip(S) of 

several experiments for random-loose sphere packings in 
two-component steel-ball mixtures have been presented. 
These are operationally defined as the random packing 
fractions obtained when pouring spheres into a container 
without subsequent densification through shaking. One 
observes that both Aip c (8) and the data follow the same 
trend. Note that we get negative values for 0.65 < 6 < 1. 
In Ref. (4l), no such effect is discussed, but it is reported 
that there seems to be no observable change in the data. 
There is no precise theoretical definition of the concept 
of random loose packing. Nevertheless, the reported val- 
ues can be taken as a quantization of a mixing effect, i.e. 
of modifications of the random cage structure. The fact 
that the variation in Atp c with 8 agrees with these ex- 
perimental findings supports the conclusion that MCT is 
able to capture the change in the average cage structure 
induced by the presence of the second component. 
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FIG. 3: Total structure factor, S(q) = ^2 a/3 S a p(q), cal- 
culated in Percus-Yevick approximation for binary mixtures 
with tp = 0.515, xb = 0.2 and three values of S. 



The results from above suggest that the change of the 
glass-transition point with composition can be under- 
stood by looking at the geometrical structure of the sys- 
tem. This information is reflected by the static structure 
factors, which comprise the relevant input for the MCT 
vertex in Eq. (Id). In particular, it is understood that 
the g-vector region around the first sharp peak in S(q) 
is important for explaining the MCT glass transition. 
Fig. U shows this region for the total structure factor, 
S(<l) = E Q ^ Saf}{q), at fixed packing fraction ip = 0.515, 
composition xb = 0.2, and different 8. One notices two 
trends caused by decreasing 8, viz. a decrease in peak 
height and an increase in its large-g wing. The interplay 
between these two trends is responsible for the shift in 
tp c . At larger 8, the increase in the wing is dominant 
and stabilizes the glass, i.e. it reduces tp c with respect 
to the one-component system. But at 8 < 0.65, the re- 
duction in peak height, equivalent to a weakening of the 
intermediate-range order, overwhelms this trend. This 
effect stabilizes the liquid, i.e. increases ip c . In all cases, 
the peak position shifts to higher q, indicating that, on 
average, particles are closer together in the mixture than 
in the one-component system; an effect typical for the 
introduction of effective attractive interactions El . 



Another way of looking at the local structure of the 
HSM is provided by the partial pair distribution func- 
tions, g a p(r). These have been obtained by numerically 
solving the Ornstein-Zernike equation in the r-domain 
using Baxter's factor function for the PY closure. Re- 
sults are shown in Fig. [|for gAj±(r) and <?ab(?")i again at 
fixed tp and Xb for various 8. Here, both quantities vary 
more or less in phase for 8 > 0.7, indicating that the lo- 
cal ordering of the one-component system is only slightly 
disturbed. One infers the B particles to be responsible 
for smaller average particle distances, thus favoring ar- 
rest of the structure. For smaller 5, the one-component 
system's structure is modified more severely, since <7aa(?") 
and (7ab(^) no longer vary in phase. Instead, "chemical 
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FIG. 4: Results within the Percus-Yevick approximation for 
the partial pair correlation functions g\\{r) and <7abM of 
binary HSM at <p = 0.516, x B = 0.2, and S = 0.9 (solid lines), 
8 = 0.8 (dot-dashed lines), and 8 = 0.6 (dashed lines). Curves 
for <7aa(?") have been shifted up by 1.0 for clarity. 



ordering" effects can be seen, and they are responsible 
for the shift of the glass transition to higher packing frac- 
tions. 

Let us stress that the variation of ip c with concentra- 
tion, while being small in total, nevertheless has a large 
impact on the dynamics close to the glass transition. This 
holds since relaxation times of the liquid in this region 
depend strongly on the distance to the critical packing 
fraction. We shall recur to this point in Sec. [y|A. 



IV. GLASS FORM FACTORS 

The spontaneous arrest of density fluctuations within 
the glass state is quantified by the glass form factors 
Fapiq)- In principle, these quantities can be measured 
in a scattering experiment via the intensity of the elastic 
line in the cross section. The diagonal elements f aa (q) 
of the normalized quantities 

fap(q) = F aP (q)/ S aa (q)S/3p(q) (10) 

have the meaning of the Debye- Waller factor for the dis- 
tribution of species a. In the limit xb — > 0, /bb(<?) is 
the spatial Fourier transform of the density distribution 
of a single localized B particle. It is the Lamb-Mofibauer 
factor /g (q) of a B particle in the hard-sphere system of 
A particles. A similar statement holds with the role of A 
and B particles interchanged, i.e., f&(q) = /aa(?, ^b — * 
1), but then the tagged particle is of the size 1/8 in units 
of the surrounding hard spheres' diameter. If the pack- 
ing fraction tp decreases towards the transition value <p c , 
the faa{q) decrease towards their critical values, faa(l)- 
These values are of particular relevance since they specify 
the so-called plateau values of the correlation functions of 
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FIG. 5: Critical glass form factors fadl) = Faa(q) I S aa (q) 
of a binary HSM with size ratio 8 = 0.8 for the large (upper 
panel) and small (lower panel) particles. The packing con- 
tributions of the small spheres are xb = (solid lines), 0.05 
(plus symbols), 0.2 (diamonds), 0.6 (circles), and 1.0 (dashed 
lines). In the upper panel, results for small q at xb — 0.01 
(crosses) and 0.1 (squares) are also shown. The dash-dotted 
line in the upper panel demonstrates the linear interpolation 
between the cases xb — and xb — 1 for qdx ^ 6. 



the liquid for states near the liquid-glass transition J42[ . 
This will be discussed further in Sec. [v| 

Fig. [| shows the critical Debye- Waller factors for small 
size disparity, 5 — 0.8, and various xb- One notices an 
increase of the values with increasing xb for almost all q. 
This result can be understood as follows. With no sec- 
ond species present, / AA (g) matches the Debye- Waller 
factor of the one-component system, f c (q), shown by 
the full line in the upper panel. On the other hand, 
as mentioned above, for xb — * 1, the quantity / AA (<7) 
crosses over to the tagged-particle quantity of a bigger 
sphere in a surrounding fluid of smaller ones, / A c (g). 
At q — ► and xb = 0, momentum conservation implies 
/aa(<? — ► 0) < 1, while for xb — * 1, particle conserva- 
tion and momentum relaxation for the tagged particle 
require / A c (<? — * 0) = 1 By interpolation, one gets 
an increase in /aa(?) with increasing xb at small q. For 
large q, on the other hand, the Debye- Waller factor in 
a one-component system is oscillating around the Lamb- 
Mofibauer factor of a tagged particle with equal diameter. 
The Lamb-Mofibauer factor in turn can be approximated 



reasonably by a Gaussian, f s (q) — exp\—(qr s ) 2 ], where 
r s is the particle's localization length ||3|. The localiza- 
tion length becomes the smaller the bigger the radius d s 
of the tagged particle is with respect to the radius d of 
the surrounding spheres fill ; in particular one gets for 
a tagged particle of diameter d s /d = 1/0.6 (1/0.8, 0.8, 
0.6) the value r c Jd = 0.041 (0.056, 0.095, 0.136). This 
implies the distribution of the / aa (<7)j gi ven hi the limit 
xb — > 1 by f s ' c (q) with d s jd = 6, to be broader than that 
in the limit £b — > 0, given by f c {q)- Therefore the width 
of the distribution / AA (q) has to increase for S < 1 as Xb 
increases from zero to unity. This is demonstrated in the 
upper panel by the dash-dotted line. It represents a sim- 
ple interpolation, f AA (q) ~ f c (q) + (f A c (q) - f c (q))x B 
for Xb = 0.6 and q > 6/d A . 

The change of /bbGz) can t> e understood along the 
same line of reasoning. But one has to notice that in 
this case the localization length of a smaller sphere in a 
surrounding of big ones matters. In particular, one has 
/ BB (<7,iB — * 0) = / B ' C ('?/'5)- This yields a width of this 
distribution smaller than the one of the /g^g, £b — ► 1) = 
f c {q/S). Such an effect can be seen in the lower panel of 
Fig. g for qdx ^ 7. The crossover is naturally given by the 
size of the A particles. Based on the above argument, one 
expects at smaller q the inverse trend. But this is only 
found in /bbO?) f° r 5 ~ ^ 7. Instead one notices that 
for all £b < 0.6, the /bb(<Z) follow closely the result for 
xb = 0, i.e. they are still close to unity at small q. This is 
a consequence of the normalization chosen here, since it is 
dominated by a change in SbbI^) at small q. It could be 
eliminated when discussing e.g. matrix-normalized quan- 
tities, f(q) = S~ 1/2 (q)F(q)S~ 1/2 (q), where the normal- 
ization properly accounts for the overall change in S(q). 

The above argument only depends on the fact that 
5 < 1, but not on the precise ratio of localization lengths. 
Thus it is quite general in binary HSM. Fig. ^| shows the 
scenario for S = 0.6, i.e. for a larger size disparity, and 
indeed one recognizes the same trends as above. Here, 
the deviations of / BB (3) from the tagged particle's / B ' c (q) 
set in faster with increasing xb than it was the case for 
S = 0.8. But one has to keep in mind that for smaller 6, 
the changes in xb induce larger changes in the number 
concentration xb, cf. Eq. (|§|). The description of / AA (<?) 
as a simple interpolation between / A ' c (?) and f c (q) as 
explained above is notably worse, indicating that this 
simple picture quantitatively only works for S not too dif- 
ferent from unity. Also more pronounced in this case are 
the changes in / BB (<7) for small q, going back to the same 
reason as outlined above. Let us note in addition that 
for both 5, the trend noticed for the diagonal elements is 
also found for / AB (q) = Fab (q) / Sab (q) , provided one is 
sufficiently far away from those q where a divergence due 
to vanishing S , A b('?) occurs. 

Macroscopic mechanic stability of the system is char- 
acterized by the elastic moduli. The liquid exhibits a 
longitudinal elastic modulus given by the structure fac- 
tor through Ml = g(kBT)J2 af3 ^S^(q -> 0)x p @. 
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FIG. 6: Same as Fig. n, but for size ratio 5 — 0.6. 



In the glass, the longitudinal modulus Ml is larger, 
Ml = M° + <5Ml, due to the arrest of the structure, 
Eq. (|^) . Figure |?] shows the results for the binary HSM 
at the transition points for 5 = 0.6, 0.7, and 0.8, and also 
the critical shear modulus, Eq. (K5|). All quantities are 
shown in units of (gh^T) in order to more clearly reveal 
the effect of composition change. Note that the total den- 
sity g of the system increases and superimposes a rise in 
the moduli one could call an "ideal mixing" contribution. 
This ideal mixing value is given by the one-component 
values, SMI ~ 56.9 and M£ w 18.3, shown through 
dashed lines in Fig. [jj. At intermediate xb , strong devia- 
tions from ideal mixing occur. For all S investigated here, 
the moduli decrease below their one-component values, 
indicating that the system becomes softer upon addition 
of smaller spheres. The effect increases with decreasing 
5 and it is of the order of 40% for 5 = 0.6. It is partly 
connected with a corresponding increase in compressibil- 
ity, k = 1/M°. Indeed, one observes for given S minima 
in all three quantities at roughly the same xb- Let us 
nevertheless point out that apart from this thermody- 
namic contribution to the softening of the glass, mode- 
coupling effects still are necessary to explain the moduli 
for S = 0.6. This can be inferred from the crossing of the 
Mf>- and SM^-versus-XB curves that is absent in M£. 
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FIG. 7: Isothermal longitudinal elastic modulus M®, increase 
of the longitudinal elastic modulus at the transition points 
<5M£, and transversal elastic modulus M-jf at the liquid-glass 
transition points in units of qIcbT as functions of the packing 
contribution xb of the smaller particles for size ratios 5 — 0.8 
(circles), 0.7 (squares), and 0.6 (diamonds). The dashed lines 
indicate the corresponding ideal mixing values evaluated for 
the one-component system. 



V. DYNAMICS 

Close to an ideal glass transition, the essential aspects 
of the dynamics are described by a universal scenario. 
This scenario has been discussed comprehensively for 
one-component systems |35| . The results of Ref. [^7| as- 
sure that these universal results are shared by the dynam- 
ics of the HSM. In particular, a two-step decay process 
arises, with plateau values given by the critical glass form 
factors, faf}(<i)i and power-law relaxations towards and 
from the plateau, governed by anomalous exponents. If 
the total packing fraction is increased towards the criti- 
cal value ip c with other parameters kept fixed, a drastic 
increase in the relaxation time r a of the slowest decay 
process is obtained that is typical for glass-forming liq- 
uids. 

In this section, we shall focus on the general, but non- 
universal features of the glassy relaxation in the binary 
HSM. To demonstrate the changes induced by different 




FIG. 8: Normalized partial density correlation functions 
4>aa{q,t) — $aa{q,t)/S aa {q) for a = A, B of a binary HSM 
with size ratio 8 = 0.8 and different packing contributions of 
the smaller particles xb for fixed packing fraction ip — 0.515. 
The wave vector is qd,A ~ 5.4. The unit of time here and in 
the following figures is chosen so that the short time diffu- 
sivity is D% = 0.01/d a . Filled diamonds mark the intersec- 
tion of the decay curves with the plateau value /° a (q) . The 
open diamonds mark a-relaxation-time scales r a (q) defined 
by 4> aa (q,T a (q))/ 7Sc(g) = 0.1. 



compositions, let us investigate a horizontal intersection 
of the transition diagram of Fig. [l] and consider a change 
of the composition xb for the total packing fraction fixed 
at ip = 0.515, i.e., at a value slightly below the glass 
transition of the one-component hard-sphere system. As 
above, the two different general scenarios, large and small 
size disparity, shall be demonstrated using the values 5 = 
0.6 and 6 — 0.8, respectively. 



A. General Features 

The dynamics for xb < 0.2 is demonstrated by Figs. |8| 
and H for the AA and BB correlation functions. We chose 
the wave vector q = 5.4/c?a below the peak in /aaC?); 
it corresponds roughly to the one used in the light- 
scattering experiment of Ref. [jjj. The normalized cor- 
relators <fiaa{q,t) = 4>aa(.q,t)/ Saa(q) cross their plateau 
values faa(q) a * certain times, say t a (q), 4>aa(q,t a (q)) = 



FIG. 9: Normalized partial density correlation functions 
4>aa{q,i) as in Fig. ^, but for the size ratio 5 = 0.6. The 
dashed lines in the upper panel show the short-time approxi- 
mation according to Eq. (|l4J) for xb = and 0.2, from left to 
right. 



faa that are marked by filled diamonds in the figures. 
Close to the transition, the correlators are close to this 
plateau for a large time interval. This is a manifesta- 
tion of the cage effect. In a leading order approximation 
for tp c — ip tending to zero, the time scale t a (q) neither 

The independence of a is 



depends on a nor on g 42 



demonstrated to a good approximation in the figures. As 
explained in connection with Figs. |^ and ||, the plateau 
increases with increasing xb and the increase is more pro- 
nounced for the larger majority particles A than for the 
smaller minority particles B. 

The decay of the correlators below the plateau is re- 
ferred to as the a process. A characteristic time scale, 
r a (q) , for this process shall be defined by specifying 90% 
of the decay: <j) aa {q, i~ a {q)) = 0.1 /J„(g), These times are 
marked by open diamonds in the figures. For 5 = 0.8, 
Fig. |8| demonstrates that the a-relaxation scale increases 
with increasing xb- This reflects the fact that with in- 
creasing xb the state corresponds to a smaller distance 
from the transition point (compare Fig. [l]). 

The scenario for 5 = 0.6, exhibited in Fig. ^, appears 
more subtle. In this case, the glass-transition diagram 
implies the distance to the transition to increase with in- 
creasing Xb, thus leading to faster decay on the a-time 
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FIG. 10: Sum fli m (g, t) of the partial density correlation func- 
tions & a a(q, t) at wave vector qd\ = 5.4, weighted according 
to Eq. (|ll|) with scattering amplitudes b a (q) as specified in 
Eq. (112). The packing fraction is kept constant at tp = 0.515, 
and ib = 0, 0.05, 0.1, 0.2 as indicated by the labels. The 
upper panel shows the results for size ratio S = 0.6, the lower 
one for S — 0.8. 



scale. But the effect of increasing plateau values with 
increasing xb was seen above to occur for all S. The 
combination of these effects leads to a crossing of corre- 
lators, as has also been observed in experiment 

In a dynamical light-scattering experiment one does 
not measure the t) directly. Rather, one measures 

a sum weighted with the scattering amplitudes b a (q) EH) , 



(11) 



a/3 



Here, Af q is some normalization constant chosen to satisfy 
cf> m (q,t = 0) = 1. It was a crucial point in Ref. 0] to be 
able to vary the b a {q) without altering the dynamics. 
Thus, three independent measurements of (j) m (q,t) could 
be used to invert Eq. (^lj) and therefore to determine the 
three distinct functions & a j3(q,t). The latter are better 
suited for a comparison with the theory. But let us also 
demonstrate the dynamics for a typical example of the 
directly measured function 4> m (q,t). If one assumes the 




FIG. 11: a-time scale ta defined through 0aa(?,ta) = 
0.1 f AA (q) f° r qd-A = 5.4 at packing fraction ip — 0.515, and 
5 = 0.8 (circles), 0.7 (squares), and 0.6 (diamonds). The lines 
are guides to the eye. 



colloidal particles to be uniform spheres, one gets [I45J 



b a (q) oc 



{qd a f 



sin(gd Q /2) - ^cos(gd Q /2) 



(12) 



Figure |Ty shows the results for 5 = 0.6 and S = 0.8 
at qd,A = 5.4. The same qualitative picture as dis- 
cussed above for the </>aa((?,£) correlator arises, yet the 
increase in plateau values is less pronounced. The reason 
is a destructive interference effect in Eq. ([n]) caused by 
®AB{q,t) < 0. This holds especially for S — 0.8, and also 
for smaller wave vectors. Nonetheless, some increase re- 
mains in all cases, and one should be able to see this in 
experiment. One could be tempted to analyze such data 
in terms of a one-component model. However, this would 
be misleading. For a one component system, the observed 
increase of the plateau f m ' c {q) = J2 a p b »(q)h(^) F ap(l) 
would imply that the system becomes stiffer upon in- 
creasing the contribution of smaller particles. But we 
have seen above from a discussion of the mechanical mod- 
uli that the opposite is the case. 

Figure |ll| exhibits a-relaxation scales ta for the larger 
particles as a function of mixing. It corroborates the 
picture suggested by the glass-transition diagram. Since 
the a relaxation close to the glass transition varies as 
r a ~ (ip c ~ <y9)~ 7 , 7 > 2.5, the variations of ta are much 
more pronounced than those of (p c . Note that the val- 
ues of ta for different S do not necessarily coincide at 
Xb = 1- MCT predicts all a-relaxation times T a (q) to be 
coupled. Thus, the qualitative picture demonstrated in 
Fig. [ll] will also hold for the a-relaxation scales of other 
experimental quantities such as the viscosities or inverse 
diffusivities. Nuclcation rates are also affected by the 
diffusivities; thus Fig. [ll] demonstrates a possible reason 
for nucleation in binary mixtures to vary strongly with 
changes of the composition. 
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B. Dynamics Close to the Plateau 

One notices in Figs. |s| and ^| a trend for the relaxation 
onto the plateau value. This part of the curve, which 
deals with the onset of structural relaxation, displays a 
slowing down of the relaxation with increasing xb for 
both cases considered for 8. In principle, the relaxation 
in this time window is a result of both structural and 
transient relaxation. In a leading approximation the lat- 
ter is given by 



=cxp[-q 2 D(q)t] 5(g) 



(13) 



with the matrix D{q) of short-time collective diffusion 
constants, D(q) — (q 2 T(q))~ 1 S(q)^ 1 . In particular, for 
a binary mixture this yields 



0AA(q,t) = l-q 2 D'(q)t + 0(t 2 ). 



(14) 



where D'{q) = x(q)D(q) with x{q) — S(q)/ 'Saa(q), an d 
D(q) and 5(g) are the diffusion constant and the struc- 
ture factor of the one-component system, respectively. 
It has already been noticed in Ref. |jj that x(g) < 1 
for small q. Thus, one expects a slowing down of the 
short-time diffusion due to mixing in the limit of small 
q. For the wave vector discussed here, however, the ef- 
fect is small: with S = 0.6, ip = 0.515, and qd>A = 5.4, 
one gets x w 0.82 (0.76, 0.78) for x B = 0.05 (0.1, 0.2). 
The approximations resulting from Eq. (14j) are shown for 
xb = and 0.2 as dashed lines in Fig. p. One infers that 
this describes the dynamics only for ^aa(<7,£) > 0.98. 
Note that x(q) is not monotonous in xb, but the men- 
tioned increase in the stretching of the short-time relax- 
ation with increasing £b is. Furthermore, at still larger 
wave vectors, one has x(q) > 1 as x(q — > oo) = X/xa, 
yielding faster short-time diffusion in the mixture. But 
the slowing down of the relaxation towards the plateau 
persists also for large g, as can be inferred from the nu- 
merical solutions. Thus we conclude that the change in 
the short-time diffusion coefficients is not sufficient to 
explain the observed effect. 

Let us now focus on the structural relaxation as de- 
fined in Sec. [n|A. Figure [j^ presents solutions of Eqs. (||) 
for 5 = 0.6 and different xb at fixed tp, together with the 
solutions reproduced from Fig. ^. The long-time parts of 
corresponding curves can be scaled on top of each other, 
since there the dynamics depends on the short-time be- 
havior only through a scaling time to, as is demonstrated 
for the xb — 0.2 curve. For other xb, the same obser- 
vation is valid. Nevertheless, we have applied the same 
rescaling as used for xb — 0.2 instead of matching to 
and t* independently for different xb- This is done in 
order to also demonstrate the drift of the scaling time 
to{xB) with xb- At short times, all structural relax- 
ation curves follow the same asymptote i -1 / 3 , and one 
notices that they deviate from one another at roughly 
t = lOt*. This demonstrates that the increase of the 
stretching in the initial decay with increasing xb , exhib- 
ited for 2.5 < log 10 (f/f*) < 5.5, is a result of structural 
relaxation rather than transient dynamics. 
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FIG. 12: Structural relaxation dynamics (solid lines) as de- 
fined by Eqs. Q for a binary HSM with S — 0.6 and ip = 
0.515, with small particle packing contributions xb = 0.2, 0.1, 
and 0.05 as indicated. The dashed lines are the solutions for 
the same parameters of the general MCT equations, Eqs. (|l|), 
with the time scaled to match the structural-relaxation solu- 
tion at long times for xb = 0.2. 



In order to achieve a deeper understanding of the con- 
clusions concerning the initial part of the structural re- 
laxation processes demonstrated above, recall that one 
can derive scaling laws for an analytical description of 
the correlators near their plateau values. This has been 
discussed comprehensively for the one-component system 
in Ref. |35| . The theory is based on the observation that 
the liquid-glass transition is described by a A% bifurcation 
of Eq. ( paj ) for the glass form factors. It is straightfor- 
ward to generalize the theory for one-component systems 
to the case of interest here Q. Let us merely note the 
basic results necessary to understand the following fig- 
ures. 

For the asymptotic expansion, one identifies a small 
parameter a and, connected to it, a time scale t a = 
toWl^ 1 ^ 2 ^- Here, the critical exponent < a < 1/2 is 
one of the nontrivial exponents of MCT that is calculated 
from the mode-coupling functional at the transition via 
the so-called exponent parameter A, A = T(l — a) 2 /T(l — 
2a). The separation parameter a is also calculated from 
the mode-coupling functional, and is a smooth function 
of the control parameters that vanishes at the transition. 
The conditions a > and a < characterize glass states 
and liquid states, respectively. Setting t — tjt a , one ob- 
tains an expansion in the small quantity 

3>(q,t)-F c (q)=H(q)^\g(i) 
+ H(q) [\a\h(t) + av] + K(q) [\a\g(t) 2 - o /(I - A)] 
+ K{q)a/ (1 - A) + 0(\a\ 3/2 ) . (15) 

Let us first explain the leading order contribution, which 
is given by the first line of Eq. (|l^). It demonstrates 
the so-called factorization theorem, in that it splits the 
wave-vector and control-parameter dependence off from 
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FIG. 13: Exponent parameters A corresponding to the points 
shown in the transition diagram, Fig. [[]; symbols indicate 8 = 
0.6 (diamonds), 0.7 (squares), and 0.8 (circles). The lines are 
guides to the eye. 
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FIG. 14: Normalized critical amplitudes /iaa(<?) = 
H A A(q)/S A A(q) for S = 0.6 and x B = 0.0, 0.1, and 0.2 as 
indicated. 



the time-dependence. It is the critical eigenvector H{q) 
introduced above that governs the former. The latter is 
given by a master function g(t) that is the solution of 



d r { 

- / g(t - t')g(t') dt' = Xg{t) 2 + sgna , 
at Jo 



(16) 



obeying g(i — > 0) ~ (t/t )~ a . The shape function g(i) 
does not depend on the details of the mode-coupling ver- 
tices, but only on the exponent parameter A. Thus, the 
factorization theorem predicts that all correlators for all 
models resulting in the same A can be rescaled to have 
the same shape, given by the master function g(t). 

The corrections to the specified scaling law are given 
by the terms of order a in Eq. (15). They consist of 
a part parallel to the critical amplitude H(q), where 
a new correction-to-scaling shape function, h(t), and a 
constant determined by the details of the mode-coupling 
vertices, v, appear. In addition, two correction am- 
plitudes, K(q) and K(q), to be evaluated from the 
mode-coupling functional, are introduced by the next- 
to-leading order. They explain that factorization holds 
with different quality for different correlators. One finds 
the K/±x(q) and -K"aa(<?) to show the same qualitative 
variation with q in the HSM considered here as in the 
one-component case discussed in Ref. The only pa- 
rameter that cannot be calculated within this approach 
is the time scale to; it is fixed by matching the long- 
time limit of the asymptotic solution at the critical point, 
& c (q,t) = F c (q) + H(q)(t/t )- a + 0(t- 2a ), to the nu- 
merical solution at long times. 

We first investigate the variation of A as a function of 
the composition, shown in Fig. [l^. The exponent pa- 
rameter is larger than the value found for the pure hard- 
sphere system, X(xb = 0) = X(xb = 1) = 0.736. It 
exhibits a maximum smaller than 0.8 for S > 0.6. As 
a result, the critical exponent decreases relative to the 



value a = 0.311 for the hard-sphere system. In partic- 
ular we get A = 0.752 (0.778), and from this a = 0.304 
(0.291) for xb = 0.1 (0.2). As a consequence, the stretch- 
ing of the decay towards the plateau increases somewhat 
with increasing xb and decreasing 6. But this effect is 
rather small and cannot explain the slowing down effect 
specified above. 

Figure [u] shows the critical amplitudes H(q) in the 
case 8 = 0.6 for the A A correlator. The normalized 
quantity /iaa(<?) = H aa (?) / S AA (?) was chosen to match 
the representation of Figs. |§| and ^|. While there is no 
general trend valid for all q, we note that at the wave 
vector q = 5.4/cJa shown above, /iaa(<z) decreases signif- 
icantly upon increasing xb- Let us emphasize that the 
region of qdA ^ 10 is the one accessible in dynamical 
light-scattering experiments on colloidal systems. Fur- 
thermore, let us add that qualitatively the same change 
with xb, although less pronounced, is observed in the 
5 = 0.8 case. The decrease of h\A yields a flattening of 
the </>(t)-versus-logi curve within the time window that 
can be described by the leading-order contribution to 
Eq. fll5j). The identified effect is further increased since 
the time scale to decreases with increasing xb- One gets 
to = 0.4408 (0.2026, 0.1385) for x B = (0.1, 0.2) and 
other microscopic parameters as given above. 

Let us turn the preceding discussion into a quantitative 
demonstration by comparing in Fig. ^ the asymptotic 
formula with the complete solution for the 4>aa correla- 
tor. The case xb — shows a typical scenario for the 
one-component system, where the first line of Eq. ( J5~5| ) 
describes over 3 decades in time of the solution, as indi- 
cated by the open diamonds. This window of the analytic 
description is expanded by the next-to-leading-order for- 
mula by about 1 decade both at short and at long times, 
as can be seen from the circle symbols. For xb = 0.1 and 
0.2, the range of validity for both the leading and the 
next-to-leading order is seen to shrink; at xb = 0.2 it is, 
including corrections, only about 2 decades. But to un- 
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FIG. 15: Asymptotic description of the normalized correlation 
functions <f>AA(q,t) for qd,A = 5.4, ip = 0.515, S = 0.6 and 
different xb as indicated. The solid lines are the full solutions 
reproduced from Fig. ^| but plotted as functions of t/to- The 
time scale t is 0.4408; 0.2026 and 0.1385 for x B = 0, 0.1, 
and 0.2, respectively. Dashed and dash-dotted lines show the 
results of Eq. (n~5[) up to order | cr | 1 / 2 and \a\, respectively. 
The diamonds (circles) mark where the asymptotic solution 
up to leading (next-to-leading) order deviates by 0.01 from 
the normalized correlator. Curves for xb = 0.1 (0.2) have 
been translated along the t-axis by 2 (4) decades for clarity. 




FIG. 16: Asymptotic description of the normalized correlation 
functions </>aa(<Z,<) for qd,A = 5.4, 8 — 0.6 and xb = 0.2 at 
different ip as indicated. Lines and symbols as in Fig. open 
symbols refer to ip = 0.515 and closed ones to <p = 0.5185. 



derstand this, one has to keep in mind that the distance 
from the critical point has increased by changing from 
xb = to ib = 0.2 for fixed total packing. Indeed, we 
get a = -0.0027 (-0.0066, -0.011) for x B = (0.1, 0.2), 
i.e. an increase in a by about a factor 4. Thus, the de- 
creasing quality of the asymptotic description is merely 
due to an increase of a. This is corroborated by Fig. jlfj , 
where the xb = 0.2 case is repeated together with a point 
closer to the transition. For xb = 0.2 and ip = 0.5185, the 
separation parameter is a = —0.0028, similar to the value 



discussed in Fig. [T| for xb — 0. Therefore, the regions of 
validity of the asymptotic expansions arc similar as well. 
Indeed, the relevant quantity specifying the range of va- 
lidity of the asymptotic expansion is not the size of the 
logarithmic time interval, but the size of the decay inter- 
val |^>aa(<7, t) — /aa(9)I- Figure s [l5| and |l6| demonstrate 
that the analytic formula, Eq. (|15[), describes the struc- 
tural relaxation in </>aa(9, t) towards the plateau /aa(<?) 
below 0.70, 0.85, and 0.90 for x B = 0, 0.1, and 0.2, re- 
spectively. This is the regime where Fig. |^ exemplifies 
the slowing down of this decay with increasing xb that 
was reported in Ref. [Q. 



VI. CONCLUSIONS 

Within the mode-coupling theory (MCT) for a binary 
hard-sphere mixture, four mixing effects have been iden- 
tified for states near the ideal liquid-glass transition. 
First, mixing suppresses intermediate-ranged ordering ef- 
fects and this leads to an increase of the small-wave- 
vector limit of the total structure factor, Fig. || Equiva- 
lently, the compression modulus of the liquid decreases. 
A similar softening of the elastic restoring forces is found 
for the moduli for compression and shear of the glass 
near the transition points, Fig. ^. Second, an apparently 
opposite phenomenon is exhibited by the increase of the 
Debye- Waller factors, i.e. a stiffening of the glass with re- 
spect to spontaneous density fluctuations, Figs. || and ||. 
This means primarily that mixing for fixed packing leads 
to better localization of the particles. The third effect is 
closely related to this, viz. a stiffening of the cages of the 
localized particles upon changes of composition. These 
changes are described by the critical amplitude, which 
decreases upon mixing, as shown in Fig. [l4|. The uni- 
versal MCT formula for the initial part of the structural 
relaxation, Eq. (filf), shows that this leads to a slowing 
down of the short-time part of the glassy dynamics, as 
demonstrated in Figs, p] a nd ^ and discussed quantita- 
tively in Figs, [l^ and gl| The above described second 
and third mixing effects have been identified originally 
in experiments on colloids 

The fourth general effect concerns the scale for the 
long-time relaxation, i.e. the scale for hydrodynamic phe- 
nomena like diffusion, or, more generally, for the a- 
relaxation processes of the liquid. Two scenarios are 
found as exhibited by the liquid-glass transition diagram, 
Fig. |l|, or by Fig. [n]. For small size disparity, mixing sta- 
bilizes the glass state. As described above for the initial 
part of the structural relaxation, also the final part of 
the decay is slowed down upon mixing. This is shown 
in Fig. H for the size ratio S = 0.8. However, for larger 
size disparities, an entropically induced plasticizing effect 
is found. Due to mixing, the glass state is destabilized 
and the a-relaxation times decrease. As a result, the 
</>(t)-versus-log(i) diagrams cross upon mixing as shown 
in Fig. |^ and observed for 6 — 0.6 in the experiments of 
Ref. m. 
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In summary, our work demonstrates that MCT can 
explain qualitatively the mixing effects on the glassy dy- 
namics of colloids observed for the size ratio 5 = 0.6. 
A quantitative comparison of the results of our theory 
with the data of Ref. jl| will be discussed in a subse- 
quent publication. Our theory suggests to also carry out 
experiments for a size ratio near 5 — 0.8 since a different 
scenario is predicted for that case. 

It can be expected that the results of our theory will 
also be of some relevance to experiments on glass-forming 
binary metal alloys. The formation of metallic glasses can 
to some extent be understood by treating the constituent 
atoms as hard spheres, which will then be all of similar 
size ]47[ . Even though in this paper we have dealt only 
with Brownian short-time dynamics relevant for colloidal 
suspensions, it is known that the long-time phenomena 
connected with the glass transition are the same for New- 
tonian dynamics ]28| , p9| . Recently, the concentration de- 
pendence of the critical temperature T c was discussed for 
a computer simulation of a Coioo-zZr-r model p8[ . This 
simulation used fine-tuned pair potentials to model the 
metallic glass former; but if one estimates the size ratio 
of a corresponding hard-sphere mixture from the atomic 
radii of Co and Zr one gets S w 0.78. Indeed, in the 
simulation the curve T c (x) was found to have a maximum 
at intermediate x, i.e. the glass transition was found to 
occur at smaller coupling strengths. But this corresponds 
to a decrease of ip c in the HSM model. A similar reason- 
ing holds for computer-simulated Ni-Zr melts [ p^[ . 

Let us add some remarks on the results derived by 
Harbola and Das E4|. Their equations, as opposed to 



the ones studied in this paper, predict for the glass-form 
factors f a p(q) and for the critical packing fraction of 
the glass transition ip c a dependence on the mass ratio 
™a/tob of the two species. This result appears surpris- 
ing because one should not expect the equilibrium re- 
sults for a classical system to depend on the particles' 
inertia parameters. It is obvious that the limit of van- 
ishing concentration .tb has to reproduce the bell-shaped 
Lamb-Mofibauer factor for the glass-form factor of the 
minority species, as discussed above in connection with 
Fig. 5. This result is not obtained in the theory of 
Ref. which, as a consequence, does not reproduce 
the experimental finding of an increase in the correla- 
tor's plateau values upon mixing [Q. Furthermore, the 
theory of Ref. predicts a much larger increase of the 
critical packing fraction (p c upon mixing than measured 
jjj. Indeed, it predicts that the glass transition can dis- 
appear completely if the size ratio 6 is smaller than a 
critical value. This result seems implausible, since there 
is no obvious mechanism which prevents the large parti- 
cles from becoming a glass upon increasing the density. 
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